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Abstract 

We study the behavior of spatially homogeneous brane-worlds close to 
the initial singularity in the presence of both local and nonlocal stresses. 
It is found that the singularity in these brane-worlds can be locally either 
isotropic or anisotropic. We then investigate the Weyl curvature conjec- 
ture, according to which some measure of the Weyl curvature is related 
to a gravitational entropy. In particular, we study the Weyl curvature 
conjecture on the brane with respect to the dimensionless ratio of the 
Weyl invariant and the Ricci square and the measure proposed by Gr0n 
and Hervik. We also argue that the Weyl curvature conjecture should be 
formulated on brane (i.e., in the four-dimensional context). 

1 Introduction 

Recently there has been great interest in brane-world cosmological models P] 
2:, particularly in an attempt to understand the dynamics of the universe at 
early times. Brane-world models have a different qualitative behavior than their 
general-relativistic counterpart f2IIJ[5JE!, especially at high energies when the 
energy density of the matter is larger than the brane tension and the behavior 
deviates significantly from the classical case. However, in spite of this interest, 
the initial singularity seems to be very little understood in these brane-world 
models. The earliest investigations of the initial singularity, which used only 
isotropic fluids as a source of matter, suggest a matter-dominated isotropic 
singularity [S1I3IH1EI - However, later work using anisotropic stresses indicate 
that the initial singularity need not necessarily be isotropic ^UJ- These issues 
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are also linked to whether or not the brane-worlds exhibit a chaotic singularity 
(see e.g. HI]). 

There are a number of important issues that need to be addressed in early 
universe cosmology. Although inflation is a possible causal mechanism for ho- 
mogenization and isotropization, there is a problem in that the initial conditions 
must be sufficiently smooth in order for inflation to subsequently take place [7| . 
This problem of initial conditions in inflation might be alleviated by a regime 
of chaotic mixing |12U13lH3| : consequently the possibility of chaotic behaviour 
in brane- world models is of interest However, more importantly, an 

isotropic singularity in brane world cosmology might provide for the necessary 
sufficiently smooth initial conditions, and might in turn be explained by entropy 
arguments and the second law of thermodynamics |17| . The study of the initial 
singularity in cosmology may consequently be closely related to the notion of 
"gravitational entropy" and the arrow of time |15I16| . A more precise statement 
of this is the Weyl Curvature Conjecture (WCC) according to which the Weyl 
tensor should be related to the entropy of the gravitational field |17lll8llT??ll2()| . 

Therefore, in this paper we will investigate the WCC for brane- world models. 
In order to do this we have to carefully study the singularity of the brane- world 
models in the presence of general fluids with anisotropic stresses (such an inves- 
tigation is lacking at present). An important conclusion from this investigation 
is, as we will show, that there exist two local past attractors for these brane- 
worlds; one isotropic past attractor and one anisotropic past attractor. For the 
isotropic singularity, the measure P defined by the ratio of the Weyl invariant 
and the Ricci square, 



decreases to zero as the singularity is approached, consistent with the WCC. 
The measure S [T^IITT?] . defined on the brane by 



behaves according to the WCC for both of these singularities. 

We note that the WCC for brane-worlds involves the five-dimensional Weyl 
tensor; additional terms arising from the projection of the five-dimensional Weyl 
tensor onto the brane should be included in the brane-world scenario. However, 
we will argue that the WCC should be formulated in the four-dimensional con- 
text. 

2 The equations of motion 

We consider spatially homogeneous brane (or branes where the spatial deriva- 
tives can be neglected). The evolution equations on the brane are as fol- 
lows |2T1I^] . The Friedmann equation, 
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the shear propagation equations, 
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Raychaudhuri's equation (9 = 3H) 
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the dark energy propagation equation, 



w + -ew + o^v^ 
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Here, if is the Hubble parameter; o"^„ is the shear tensor; ^TZ^ is the 3- 
curvature; p the energy-density; p the isotropic pressure; 7r M1/ the anisotropic 
stress tensor; U and Vp V are the nonlocal dark energy and the bulk graviton 
stress tensor, respectively. The nonlocal energy-flux Q M has to vanish due to 
the shear constraint equation. Note that there are no propagation equations for 
the nonlocal bulk graviton stress tensor, V^. 

For a FRW brane near the initial singularity, the Friedmann equation gives 



H 2 cx p 2 , 



(7) 



and thus H, p cx t -1 . The three-volume, v, goes as v cx t 1 ^. We will primar- 
ily be interested in spatially homogeneous and spatially flat branes. The FRW 
universe corresponds to an equilibrium point of the equations of motion when 
written as a dynamical system, and hence may serve as a local past attrac- 
tor. The behavior of the brane-worlds close to these equilibrium points, which 
have the important property that H cx t , are reduced to linear algebra by 
linearizing the set of equations with respect to the equilibrium point. 

In the following we will assume that the brane-world contain two types of 
fluids; an isotropic fluid obeying the barotropic equation of state pi = (7— l)pi, 
and a radiation type of fluid causing the anisotropic stresses ir^. The fluids 
separately obey the energy-momentum conservation equations 



pi + 76p l 
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3 Isotropic or anisotropic initial singularity? 

In an earlier paper ^U] it was shown that the initial singularity for a magnetic 
brane-world could be either locally isotropic or anisotropic. The question we ad- 
dress here is: What happens in the case with more general types of stresses? Will 
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the same behavior persist? Will there exist one anisotropic and one isotropic 
past attractor for more general brane-worlds? 

In contrast to the general-relativistic models, the brane-worlds will in gen- 
eral be matter-dominated as we approach the initial singularity. The question 
whether the singularity is anisotropic or isotropic reduces to the question of 
which fluid dominates at the initial singularity. If the anisotropic fluid dominates 
initially the fluid will typically drive the initial singularity to an anisotropic 
state. If, on the other hand, the isotropic fluid dominates, the universe will be 
driven to an isotropic FRW state. 

3.1 Isotropic singularity 

Let us first assume that the initial singularity is isotropic; hence we have a FRW 
singularity with H, pi oc t _1 so that the three- volume expands as a power-law 
v oc i 1 / 7 . Moreover, isotropy implies that the anisotropic stresses and fluid have 
to be small in the sense that ir^, a^, p r <C pi. 

To check the stability of the isotropic singularity we first investigate the 
stability of the ratio of p r /pi- For the isotropic singularity to be stable for this 
combination of fluids, we must have p r /pi — > as t — > 0. Using eq. JHJ, we get 
the following requirement for the isotropic singularity to be past stable (using 
also that 7r M „, cr^u, p r < p t ) 

1 > 4/3. (9) 

It is quite easy to understand this requirement intuitively. Close to isotropy 
the anisotropic fluid behaves as isotropic radiation. If the isotropic fluid has 
7 < 4/3, then the anisotropic fluid will be dominant initially. This would, 
however, drive the universe into an anisotropic state and hence the isotropic 
FRW would be unstable. If, on the other hand, the isotropic fluid is stiffcr 
than radiation, then the isotropic fluid dominates initially and the isotropic 
singularity is stable. Close to the isotropic singularity the shear will behave as 
a 7 = I fluid and will therefore also be suppressed near the initial singularity. 
We note that 7 > 4/3 for scalar field matter close to the singularity (which acts 
effectively as stiff matter with j e ff = 2). 

The only thing we have assumed in the above heuristic analysis is that the 
dominant term in the Friedmann equation is the pf-term. The dark-energy 
U behaves close to isotropy as classical radiation and will therefore be sub- 
dominant to the anisotropic radiation fluid (which come with quadratic terms 
in the Friedmann equation) . The nonlocal stresses V^ u have the same source as 
the dark-energy (they both are different components of the Weyl tensor in the 
bulk) and hence one should expect that the nonlocal stresses are sub-dominant 
as well. The above stability analysis is therefore quite general and encompasses 
many different types of brane- world scenarios. 

To find the asymptotic behavior of the shear we have to specialize to a certain 
type of anisotropic fluid. We assume that the anisotropic stresses are of the form 
= C^ v p r , where is a matrix of constants (or slowly varying functions). 
This choice includes the magnetic field as a special case and may give us a useful 
indication of the asymptotic behavior of the shear in the presence of stresses. 

Consider the isotropic singularity with 7 > 4/3. Assuming that behaves 
approximately as the dark-energy term we can neglect the nonlocal stresses on 
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the right-hand-side of the shear propagation equations. To lowest order the 
equations are 

<j(nv) + 9°> oc pip r CX < 3^. (10) 
This equation can be integrated to give to lowest order 

= O (i - ^) . (11) 
3.2 Anisotropic singularity 

This case is more difficult than the isotropic one since the anisotropic singular- 
ity is almost entirely determined by the anisotropic matter. Also, the nonlocal 
stresses come into play and are difficult to deal with in full generality. The 
main reason for this is that from a brane point of view, there are no evolu- 
tion equations for V^v In principle, V^v should arise as a solution to the full 
five-dimensional theory, but we shall assume V^v = in the specific example 
(equilibrium point) discussed below. However, by a continuity argument we can 
expand the results to V^ v = lAD^, where D llv D tlv is sufficiently small. 

We start out with the assumption that there exists an anisotropic equilibrium 
point in the past. We can then give a criterion that must be fulfilled if the 
equilibrium point is past stable. For an important class in which there exists 
such an equilibrium point, we shall show that the equilibrium point is past stable 
if and only if this criterion is satisfied. 

We assume that there exists an anisotropic equilibrium point such that 
ti>„, (J^v fx t , and pi <C p r - We also assume the asymptotic values of p r 
and are p r = r/t and = 6o/t, respectively, where r and 9q are constants. 
This implies that we have the asymptotic value 

a 1^ w constant = D (12) 

Pr 

Using eq. (JSJ we get the following asymptotic relation 

-l + ^9 + Dr = 0. (13) 

Consider the ratio R = pi/p r - For the equilibrium point to be past stable R > 0, 
so using eqs. (jHJ) and lfT3)l. we obtain 

#07 < 1. (14) 

The constant 9q is related to the volume expansion rate as v oc t e °, so eq. I|14(l 
gives an upper limit on the exponent on the power-law behavior of the volume 
close to the equilibrium point. Note that if do < 1/2, then this inequality is 
automatically satisfied for all 7 6 [0, 2]. 

So far we have only analysed the stability in the R direction. In order 
for the equilibrium point to be past stable, it has to be stable in the remaining 
variables as well. Let us now provide an example which is past stable in the other 
variables. The example is a two-parameter set of different matter configurations 
which interpolates between the magnetic brane-worlds and a brane-world with 
two isotropic fluids. Only ea. (|14f> then remains to be satisfied. 
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Given an anisotropic stress in the diagonal form 

Tr^iy = Prdiag (o, x + %/%, x - V3y, —2x\ , (15) 

where x, y are constants and we assume that the nonlocal stresses vanish, V^v = 
0. We will also assume that the dominant energy condition 

\Pj\ = Y + 7r " ~ Pr ^ 

is satisfied for the fluid p r . This implies that (x,y) are restricted to a compact 
subset of M 2 . 

We define the following constants 

A = 2(x 2 + y 2 ) +x 3 - 3xy 2 

S 2 = (2x + x 2 - y 2 ) 2 + (2y - 2xy) 2 

L 2 = 2 [2 + 9(x 2 + y 2 )] A 

a = 36A 2 + 8£ 2 + 6L 2 



B = ^(a + L 2 ) 2 -8a(Z 2 +L 2 )-L 2 . (17) 
Then there exist a past equilibrium point given by the asymptotic values 
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3(o- B) 1 
4a i' 
(2x + x 2 - y 2 )B 1 
6Aa t' 
(2y - 2xy)B 1 



(18) 



6Aa i ' 
k 2 2 _ B{B + a) 1 
T Pr ~ 12 Aa 2 i 2 ' 

p 4 = 0, 

2U _ 9A 2 (a - Bf - 45 2 E 2 - 2Aff (5 + a) [l - 9Qr 2 + y 2 )] 1 
~ 144A 2 a 2 t 2 

Here are cr± related to a^ v via 

= diag [O, a+ + V3a-,a + - \/%(T-, — 2a + ^j . (19) 

It can be shown that these asymptotic solutions are past attractors within the 
subset V^v = 0, provided that inequality (jl4|) holds. Further, it can be shown 
that 9q < 3/4 for nonzero x,y. 

Using a continuity argument 1 we can show that there also exist a past attrac- 
tor for models with nonlocal anisotropic stresses of type V^u = UD^ V where 
is sufficiently small. Hence, there is a class of models with V^v ^ which have 
an anisotropic past attractor. How large this class is, and if this anisotropic 
past attractor exists for generic brane-worlds, needs further work. 2 

1 The argument uses the fact that finding the equilibrium points of a dynamical system is 
basically a question of finding solutions to an algebraic set of equations. Also, the stability 
analysis for these equilibrium points depends on the eigenvalues of a set of linear equations 
which are continuous functions of the C M „ and D^ v . 

2 The paper 1231 considers brane dynamics with nonlocal stresses. However, the brane- 
world does not contain an additional anisotropic fluid and thus there are no anisotropic past 
attractors. 
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4 The Weyl Curvature Conjecture 



We have seen how different brane- world scenarios have two different past attrac- 
tors. Here we will explore both of these possibilities further and check whether 
the brane- worlds behave according to the WCC (see |18lll9j for a review). 

For the isotropic singularity we need to know how the Weyl tensor approaches 
the asymptotic values. In the previous section we found that er M „ oc £~ 4 / 3 t (recall 
that 7 > 4/3 in this case). The dominant term in the Weyl curvature invariant 
will therefore be (<7 M „) 2 . Hence, 

(CtfT a C a/ s 7 «)j = t- 4 ■ O (^ (37 ~ 4) ) . (20) 

The Ricci square behaves as t~ 4 , thus 

Pj = o(tW*-*A, 

Si = O (t^- l) \ . (21) 

Hence, the isotropic singularity behaves according to both the P-version and 
the iS-version of the WCC in this case. 

For the anisotropic singularity, the behavior is easier to establish. The Weyl 
curvature behaves according to 

(C a ^ s C aPl& ) A =t- i . (22) 

The Ricci tensor evolves similarly, so 3 

Pa — constant, 

S A oc t e °. (23) 

Recall that 9o obeys the bound #o7 < 1) an d the explicit value of 9o depends 
very much on the matter content. In all of the cases we have explored here, 

e < 3/4. 

We can conclude that the measure S behaves according to the WCC for both 
the anisotropic and the isotropic singularity. 



5 A four-dimensional or five-dimensional WCC? 

So far we have considered only the four-dimensional version of the WCC, seen 
from the brane point of view. However, a question that arises for these brane- 
worlds is whether one should consider the WCC from a hve-dimensional per- 
spective or a four-dimensional perspective. Even a four-dimensional measure 
of "gravitational entropy" should have a part arising from the projection of 
the five-dimensional Weyl tensor onto the brane. Hence the gravitational en- 
tropy on the brane will have one part corresponding to the four-dimensional 
Weyl tensor, and a part from the remaining components of the projection of the 
five-dimensional Weyl tensor. These components are 

= ^CABCDn A n c g B ^g D v , (24) 
jW = i5) C AB cDn D g c a g A gl, (25) 

3 Although Pa is asymptotically a constant, its derivative can be either sign asymptotically. 
4 It should be noted that all these entities should be defined in the bulk in the limit close 
to the brane. The Weyl tensor itself experiences a delta function singularity at the brane 1241 . 
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where n A is a unit vector orthogonal to the brane, and g A ^ is the projection 
tensor onto the brane (see [SSIEHI) • These nonlocal terms give rise to "dark 
energies" and stresses that do not come from standard model particles. More 
explicitly, the tensor can be decomposed as 



£ =__L 



(26) 



and thus is, up to a constant factor, formally equivalent to the energy-momentum 
tensor of a radiation fluid with energy-density U, anisotropic stress V^, and 
energy-flux Q^. The entropy of particles is calculable from standard thermo- 
dynamics, and hence there should also correspond an entropy due to the dark 
energies. These terms are purely gravitational effects but behaves very much 
like ordinary matter when it comes to the evolution of four-dimensional space- 
time. Actually, the notion of gravitational entropy might even be more easily 
understandable for these terms than for the more mysterious Weyl tensor. 
The five-dimensional Weyl curvature invariant can be written as 

( 5)c abcd Cabcd = WcxPrtCafrs + GS^S^ + AB^ a B^ a . (27) 

The tensor B^ va can be expressed in terms of the exterior curvature of the 
hypersurfaces as [26] 

B^a = 2W {ll K v]a + I (V/jJf^ - V^K) g v]a (28) 

where the covariant derivative, V M , is the covariant derivative associated with 
the metric on the brane, g a f}- Further, from the junction conditions and the 
1*2 symmetry, we can relate the exterior curvature close to the brane and the 
energy-momentum tensor on the brane, 

K tu = ~ K »v = -\i4hy.u - \f \ (t^ - \g^T^j . (29) 

In our case, where the brane is spatially homogeneous, the components B^a 
can be replaced with time-derivatives of the matter on the brane. 
There are different versions of the WCC that can be investigated. 



5.1 The WCC: Modified four-dimensional version 

First, we will consider the four-dimensional version where we take into account 
the extra terms from the projection of the five-dimensional Weyl tensor. Only 
the four-dimensional versions of the tensors shall be used. 

For the isotropic singularity, we will consider a FRW brane for which exact 
five-dimensional solutions are known. In this case, 

(B^ a B^ a ) FRW = 0, (30) 

using the junction conditions and energy-momentum conservation of the fluid 
on the brane. The term S^S^ is, in general, given by 
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Hence, for the FRW brane 5 , 

{ 4) c aM8 Ca ^ s + Qg^ gfiv + 4 B ^ a B^ a oc U 2 . (32) 

For 7 > 2/3 the Ricci square on the brane diverges as p 2 , so 

Pfrw « t^- 2 ^ (33) 
Sfrw oc t^C 6 ?- 1 ). (34) 

For the anisotropic singularity, we get a similar behavior as in the pure brane 
case: 



^ )c abcd Cabcd 



cx t- 4 . 

A 

Pa = constant, 

S A oc t e °. (35) 



5.2 The WCC: The five-dimensional version 

The five-dimensional curvature tensors experience (^-function singularities at 
the location of the brane and thus care is needed in this case. Let us therefore 
consider the case where the volume under consideration is entirely in the bulk 
but very close to the brane. 

The measure = ^C abcd C A bcd/ (5) R MN Rmn can be used to inves- 
tigate the five-dimensional WCC in the bulk. The bulk, which is taken to be 
an Einstein space with negative curvature 6 , has ^Rab oc Ab^Qab- On the 
brane, the five-dimensional Ricci tensor also experiences a ^-function singular- 
ity. Hence, it is of crucial importance where one wants to consider the WCC. 
If we consider ' 5 'P in the bulk, close to the brane, we can use the known exact 
solution for the FRW brane. 

The FRW-brane can be embedded in an ambient Schwarzschild-AdS space 

m 

ds 2 = -f(R)dT 2 + j^- + R 2 1%] dx l dx\ f(R) = 1 - J - ^rR\ (36) 

where 7^ is the 3-dimensional spatial metric encompassing the fifth dimension. 
For simplicity, we have also assumed the closed FRW model. The induced metric 
on the brane is 



ds 2 = -dr 2 + R{rf 



^ + r 2 (d0 2 + sin 2 9d(f> 2 



(37) 



where r is the proper time of the brane, and R(t) describes the motion of the 
brane; i.e. R(t) can be identified as the scale factor a(t). We can then easily 
find the Weyl curvature invariant: 

( { ,) c ABCD Cabcd \ Q(i? -8_ (38) 
V / FRW,bulk 



5 For an anisotropic brane with shear degrees of freedom, in addition to this mode we get the 
usual shear modes from the four-dimensional Weyl invariant, ^C a ^ lS C a p 1 s ) found earlier. 

6 Recently a scalar field in the bulk has also been included, but these spacetimes have only 
been investigated at a perturbative level. 
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Hence, 



^P) oc R-\ 

FRW.bulk 



oc R- 1 . (39) 

FRW.bulk 



Thus, both measures diverges as we approach the initial singularity. 

For the anisotropic singularity, unfortunately there are not many exact so- 
lutions known to study the possible behavior. However, there are exact Kasner 
solutions HE1> and we can use these solutions to obtain 



-►oo, as t -> 0. (40) 

A,bulk 



Therefore we conclude that the measures investigated here diverge badly in the 
bulk as we approach the initial singularity. Thus the WCC as formulated above 
is not valid in the bulk. 

This is perhaps not very surprising since the physics in the bulk is very 
different than the physics on the brane; indeed, the matter in the bulk need 
not even obey the weak energy condition. Moreover, matter is confined to the 
four-dimensional brane, and the physical energy conditions are only satisfied on 
the brane. Therefore, we conclude that the WCC should be formulated in the 
four- dimensional context. 



6 Conclusion 

In this paper we have presented a thorough investigation of the initial singularity 
in brane-world cosmological models. It was shown that for a class of spatially 
homogeneous brane-worlds with anisotropic stresses, both local and nonlocal, 
the brane-worlds could have either an isotropic singularity or an anisotropic 
singularity. 

We also discussed the Weyl Curvature Conjecture, according to which the 
gravitational entropy is related to the Weyl curvature, in the context of the 
newly proposed brane-worlds. We found that the isotropic equilibrium point is 
typical in the brane-world scenario. 

For both the isotropic and the anisotropic singularity the Weyl curvature 
invariant, C af3 ~< s C a ^ 5 , the Weyl to Ricci ratio, P 2 = G a ^ 5 C a ^ s / 'R^ R^, and 
the measure S = VhP were investigated. The measure S behaved according 
to the WCC in both cases, while P showed generic increasing behavior only for 
the isotropic singularity. Hence P and S can be distinguished as measures of 
gravitational entropy. 

However, for the brane-worlds a "gravitational entropy" may also have ad- 
ditional terms arising from the higher-dimensional bulk. The nonlocal terms, U 
and Vfiv, are projections of the five-dimensional Weyl tensor and behaves very 
much like a fluid upon the cosmological evolution of the universe. Hence, these 
terms are, in some sense, a direct manifestation of the Weyl tensor in terms of 
dark energies and stresses. These terms should therefore be incorporated into 
a gravitational entropy, one way or another. However, it seems as if the WCC 
can only be formulated on the brane where it can be seen upon as the "shadow" 
of the five-dimensional theory. 



10 



Finally, we remark that brane-worlds may serve as an interesting arena for 
studying the WCC. Actually, the brane-worlds may serve as a scenario were 
gravitational entropy and its relation to the Weyl tensor may be more easily 
understood. 
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